Coulomb blockade in two island systems with highly conductive junctions 
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We report measurements on single-electron pumps, consisting of two metallic islands formed by 
three tunnel junctions in series. We focus on the linear-response conductance as a function of gate 
voltage and temperature of three samples with varying system parameters. In all cases, strong 
quantum fluctuation phenomena are observed by a log (JzbT / (2E co )) reduction of the maximal 
conductance, where E co measures the coupling strength between the islands. The samples display 
a rich phenomenology, culminating in a non-monotonic behavior of the maximal conductance as a 
function of temperature. 

PACS numbers: 73.23.Hk, 85.35.Gv, 73.40.Gk 
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The transport properties of single electron devices are 
well described by the so called sequential tunneling model 
as long as the conductances of the underlying tunneling 
contacts are small compared to the conductance quan- 
tum G K =e 2 /h 00. In most practical cases, however, 
one has to include higher-order corrections for a quan- 
titative description due to quantum fluctuations of the 
charge number on the single electron islands. This has 
been demonstrated clearly for the single electron tran- 
sistor (SET) H H II II Excellent agreement be- 
tween experimental investigations and theoretical studies 
has been established at low conductances with the aid 
of perturbation expansion including the correct descrip- 
tion of a logarithmic contribution to the linear-response 
conductance of the multichannel-Kondo type. At higher 
conductances however, perturbation theory breaks down. 
Nevertheless, by applying Monte-Carlo methods, the ex- 
perimental results can be described with amazingly high 
accuracy 0, 0, The findings on the SET indicate 
that a close match between the model Hamiltonian and 
its experimental realization exists for single electron de- 
vices. It is not clear a priori whether the good agreement 
between perturbation expansion and experiment as well 
as Quantum-Monte-Carlo numerics found for the SET 
survives if more complex arrangements of single electron 
islands are investigated. E. g. the low temperature con- 
ductance of the SET involves finite occupation of only 
two states and the mapping onto the multichannel-Kondo 
model is based on identifying those two states with a 
pseudo spin. This procedure has no simple analogy in 
general arrangements. It is worthwhile to check the range 
of validity of perturbation theory in more general cases. 
Two island systems, readily accessible both by experi- 
ment and theory, serve as a good starting point. 

In this Letter we present an experimental study of 
the linear-response conductance of two single electron is- 
lands in series, an arrangement nicknamed single electron 
pump (SEP) 0, E3, E3 and sketched in Fig. □ The 
linear-response conductance varies with the gate volt- 
ages and is bound between temperature dependent val- 



1/, o- 



V, o- 



u, 



n 



r 



-OV, 



-oVo 



FIG. 1: Schematic representation of the single-electron 
pump. The three tunneling contacts Ti, T m , and T r are laid 
out in a row forming two islands (each contact i £ {I, m, r} 
is characterized by its capacitance d and conductance Gi = 
giGx)- ni and n r count the number of electrons by which the 
left and right island charge differs from neutrality. The ar- 
rangement is biased by the voltage difference V r — Vi. The elec- 
trostatic potentials on the islands (Ui and U r ) can be tuned by 
Vi and V2 which couple directly to the islands by the capaci- 
tances C11 and CV- C21 and C\ r represent the experimentally 
unavoidable stray capacitances. 



ues Gmin(T') and G max (T). For samples in an interest- 
ing and accessible parameter range, the maximal conduc- 
tance G max obeys a non-monotonic behavior as a func- 
tion of temperature, which is reported here for the first 
time. This behavior, astonishing at first sight, is natu- 
rally explained by the sequential tunneling model in this 
Letter, giving evidence that non trivial behavior is to be 
expected in complex arrangements of single-electron is- 
lands. Beyond the validity range of the sequential model, 
we find a logarithmic correction of the conductance at 
low temperature due to quantum fluctuations. Our high 
quality data make a sensitive test of perturbation theory 
as developed e. g. in Ref. feasible. This however, re- 
quires an elaborate calculation along the line of formulae 
in Ref. 0, which is out of the scope of the work presented 
here. 

The linear-response conductance of the SEP can be 
modeled in various coordinate systems spanning the 
(I/11V2) plane. In this letter we choose dimensionless 
coordinates n x and n y defined such that the charging 



energy E ch = E cx (n x - 



- E m (n y + An + nn s ) . Here 
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FIG. 2: a) Stability diagram for the SEP sketched in Fig. [T] 
The hexagonal cells mark regions where the indicated charge 
state (m, n r ) possesses the lowest energy E c h(ni,n r ). b) Mea- 
surement of the linear-response conductance of sample 2 as a 
function of both gate voltages Vi and V2 at 150 mK. The out- 
ermost contour line indicates a conductance of 0.01 uS. The 
following lines range from 0.05 uS to 0.4 uS with a 0.05 uS 
spacing. The stability diagram (thick lines) is deformed in 
the coordinates of this figure. 

E cx = e 2 /(2C s ) with C s being the sum of all external 
capacitances in Fig.n(i-e. all capacitances except C m ), 
E m = (e 2 /2)C s /{C m C s + C s iC sr ) with C s i and C sr de- 
noting the sum of the external capacitances on the left 
and the right island separately, and k — (C sr — C s i)/C s . 
The coordinate n x is associated with the change of the 
total charge number n s = ni + n r , while n y redistributes 
the charge between both islands (Art = ni — n r ). 

Each point in the plane spanned by n x and n y can 
be mapped onto a charge ground state (nj,n r ) which 
gives the lowest possible electrostatic charging energy 
E c h(rii, n r ). This procedure divides the (n x ,n v ) plane 
into the grid of hexagonal cells depicted in Fig. |2Ji. The 
length E co / E cx of the horizontal cell boundaries is a mea- 
sure of the coupling strength, where E co = E cx +E m {n 2 — 
l)/4 cx C m . Due to the periodicity it suffices to study the 
linear-response conductance in a small exemplary portion 
of the (n x ,n y ) plane. Within the sequential model, the 
linear-response conductance of the system vanishes expo- 
nentially at low temperatures except close to the triple 
points in the (n x ,n y ) plane, where the ground-state en- 
ergy of all three adjacent states is degenerate. The con- 
ductance peaks near these points. It is worthwhile to 
mention a peculiar behavior of the SEP: Exactly at the 
triple points the low temperature conductance is constant 
and given by G /3 where Gq 1 = G^ 1 + Gy 1 + G m x . Al- 
though for T <C E co /kB the maximal conductance G max 
is temperature independent as well, the system assumes 
Gmax at a slightly different position in the n x direction 
[12j. Taking only three states into account (the occupa- 
tion probability of all other states is exponentially small 
for T < E co /k B ) we get 

G s ( g s e /3Ai5 -l\ _1 
G(n x )\ ny= _ K = - - e ° 0AE \ -L + ^ E j , (1) 



FIG. 3: Scanning electron microscope picture of sample 2. 
The two T-shaped structures in the middle are the islands. 
The inner tunneling contact is visible in the center of the 
picture. The outer contacts are turned by 90° with respect 
to the former. The gate lines are visible at top and bottom 
pointing up and down towards the inner contact. In the final 
measurement the top and bottom leads are biased in parallel. 

with (3 = l/(k B T), g s = gig r /(gi + g r ) = G S /G K , and 
AE = 2E cx (n x — 1) +E co . Depending on the ratio g s j g m 
the peak position deviates from the triple point at AE — 
for g m ^ 2g s . 

In a measurement as a function of n x and n y the con- 
ductance displays a periodic grid of peaks (grouped as 
pairs), marking the endpoints of the horizontal bound- 
aries in Fig. With rising temperature the peaks 
broaden and shift towards the center of the horizontal 
boundaries. Fig. 03 gives an example at T = 150 mK~ 
0.09E cx /kB, where thermal broadening already is effec- 
tive in merging the two separate peaks. Finally the two 
peaks merge completely and the conductance takes its 
maximal value at the mid-points, e.g. at n x = 1 and 

Tly = —K. 

We study three samples in two different layouts. 
All samples have been produced by standard shadow- 
evaporation technique from aluminum with aluminum- 
oxide barriers. The barriers for the middle and outer 
contacts are fabricated in different oxidation steps, mak- 
ing different barrier thicknesses for internal and external 
contacts possible. Sample 1 has a simple layout (not 
shown) with three contacts in a row, which is the most 
natural arrangement to fabricate a SEP. In this layout the 
serial conductance is accessible, but the conductance of 
the individual contacts remains unknown. In the slightly 
more complex structure of sample 2 and 3 (see Fig. 0J) 
each island is connected to independent leads via two 
contacts. This permits measuring the conductance of 
different serial combinations of the contacts, and thus 
the individual contact conductances can be determined. 
In the final experiments both external contacts of each 
island are operated in parallel (connected to the same 
voltage source) , and the two contacts then act exactly as 
a single contact. We could not detect any degradation 
(considering noise performance or sensitivity to external 
disturbances) for the latter samples compared to sample 
1. 

From the positions of the conductance peaks in the 
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TABLE I: Parameters of the three samples. 



samp. 


(K) 


■t^m ^co -t^c 

ks ks ks 
(K) (K) (K) 


K 


gi 9r g m 


Go 
(pS) 


1 


2.8 


5.8 1.3 5.6 


-0.018 


0.44 0.44 0.04 


1.42 


2 


1.6 


3.0 0.9 2.5 


0.10 


0.52 0.83 1.32 


10.0 


3 


1.5 


4.7 0.3 4.5 


0.0013 


0.73 0.57 0.03 


0.95 



(n x , n y ) plane at low temperatures (see Fig. [SJl we obtain 
the parameters ft and E co /E cx . In addition we measure 
the conductance in the high temperature region where 
it does not depend on either n x nor on n y . In close 
analogy to the high temperature expansion for the SET, 
it behaves as G(T) w G ((1 - E c /(3k B T)) . It can be 
shown [l^llij that the relation 
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holds, where e± = ((E co /E cx ) (ft ± 1) =p 2) / (ft =F 1) and 
E c is an experimentally determined fitting parameter. 
Tab. [I] gives all relevant sample parameters. 

To simplify the analysis of our data and facilitate the 
comparison with theoretical considerations we focus on 
the temperature dependence of three special values of the 
conductance in the (n x ,n y ) plane, namely G min , G max , 
and G m , the latter defined as the conductance at the 
center of the horizontal boundaries in Fig. (i. e. at 
n x = l,n y = —ft). In Fig. 0] we display our main find- 
ings. In addition to our experimental results we present 
the outcome of a calculation in the framework of the se- 
quential tunneling model using the parameters from Tab. 
U The techniques for such a calculation are well docu- 
mented in the literature (e.g. Q) so we do not comment 
on this calculation in detail. It requires the solution of a 
master equation which one may set up using golden rule 
rates for the inelastic tunneling events. For T -C E c /ks 
at most four (ni,n r ) states are occupied with reasonable 
probability. Restricting the master equation to these four 
states allows for an analytical solution. At higher tem- 
peratures numerical relaxation methods are used. 

The sequential tunneling model gives good agreement 
with our experimental data with the exception of G max 
at low temperature. The latter deviations are discussed 
at the end of the letter. The overall behavior (see Fig. 0} 
is governed by two scales (E co /kB and Ec/ks) at which 
G m (T) and Gmin(T'), respectively, start to increase and 
finally merge with G max (T). To get finite conductance 
through the SEP at least three charge states have to 
be occupied in thermal equilibrium. At T <C E co /kB 
this is only possible near the triple points where three 
adjacent states are occupied. Here e. g. the sequence 
(0,0) -> (1,0) -» (0,1) -> (0,0), corresponding to a 
charge transfer from left to right, occurs with finite prob- 
ability. The inverse process is equiprobable, but under 
voltage bias a net current occurs. At n x = 1 
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FIG. 4: Linear response conductance of the three samples 
as a function of temperature. Shown are G m m (o), G max (•), 
and G m (0, the conductance at n x — 0.5). The dashed lines 
are the result of an analytical solution of a 4 state model (see 
text). As solid lines the outcome of the sequential model with 
parameters from Tab.[I]is shown. For the thin solid line a term 
of the form alog(fcsT '/ 'E co ) has been added to the maximal 
conductance as calculated from the sequential model. Here a 
is a fitting parameter. 



where G m is measured, the charging energy of states 
(0,0) and (1,1) lie AE = E co above the two fold de- 
generate ground state ((0, 1), (1,0)). As a result the lin- 
ear response conductance at T <C E co /kB, as calculated 
from the sequential model, is exponentially suppressed 
since no charge transfer is possible using only two states. 
At T ~ E co /kB the states (0,0) and (1, 1) are thermally 
occupied with increasing probability leading first to an 
exponential increase of G m and finally to the merging of 
G m and G max . For n x — n y — states besides (0, 0) are 
occupied for T > E c /ks only. Thus G mm is exponen- 
tially small for T < E c /k B - 

The most striking feature of our measurements is the 
non-monotonic dependence of G max on the temperature 
found for sample 3. It is clear from Fig. 0] that the 
phenomenology is correctly described by the sequential 
model. To get more insight into the nature of the drop 
of G max at T > E co /kB we analyze G m (T), which coin- 
cides with G max in the relevant temperature range. The 
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FIG. 5: Position of the conductance maxima as a function 
of temperature. •: sample 1, ♦: sample 2, ■: sample 3. 
Lines: Predictions of the full sequential model (solid) and Eq. 
Q (three-state approximation, dash dotted) . For sample 1 the 
peak position has been used to adjust g s /g m = 5. Assuming 
9i — 9r — g-m, gives the thin lines failing to fit the data. 

four-state approximation mentioned above yields an an- 
alytical solution: 

r (4) = G s fiE co ( 2g s j3E co \ 1 

2 smh{(3E co ) \ gml-ePz-J ' 1 1 

We have drawn this function for our sample parameters 
in Fig. 21 as dashed lines. It has a distinct maximum for 
g m < g s , the position and strength of which depends on 
the ratio g s /gm- The above approximation breaks down 
at T ~ E c /kB where more than four states are occupied 
significantly resulting in a rapid rise of the conductance. 
For an experimental observation of a local minimum the 
relation E co < E c has to be fulfilled in addition to g m < 
9 s- 

In Fig. [SJwe plot the n x coordinate of the conductance 
maximum n max as a function of temperature. Again, rea- 
sonable agreement with the sequential model is found. 
For T — > the position of the maximum approaches 
the location of the triple point. At low temperature 
the conductance peak is described by Eq. ^ The maxi- 
mum can shift in either direction depending on the ratio 
9s I 1 gm- This can be used to determine g s /g m for sam- 
ple 1 {g s / 9m = 5) where this ratio can not be measured 
directly due to the simple layout (see above). 

The sequential model fails to predict the temperature 
dependence of G max at low temperature. This deviation 
can be described phenomenologically by adding a term of 
the form AG qf = a log(k B T/(2E co )). The thin solid lines 
in Fig. 01 display the outcome of a fitting procedure in a 
that minimizes the mean square deviation between the 
measured values of G max and G scq + AG q f(cv), showing 
very good agreement. We find ot\ — 45 nS, 012 = 200 nS, 
and CV3 = 33 nS for sample 1, 2, and 3, respectively. For 
the SET much the same behavior is observed and at- 
tributed to quantum fluctuations of the charge states. 
Pohjola et al. .12| analyzed the linear response of the 
SEP by renormalization group methods. They also found 
a logarithmic behavior of the low temperature conduc- 



tance, in qualitative agreement with our experimental 
result. However, quantitative results were obtained for 
certain limiting cases of special interest only. For a de- 
tailed comparison with our experiment a calculation us- 
ing our sample parameters is highly desirable. 

In summary we have presented an experimental study 
of the linear response conductance of the SEP in a regime 
where quantum fluctuations of the charge eigenstates can 
not be ignored. Depending on the ratio g s / 1 g m the SEP 
shows a remarkably rich phenomenology even within the 
framework of lowest-order perturbation theory, the so- 
called sequential tunneling model. Most strikingly, in the 
easily accessible regime E co < E cx and g m <C 2g s a pro- 
nounced non-monotonic temperature dependence of the 
conductance has been observed. Phenomena of this kind 
are to be expected in all single electron devices which 
are more complex than the SET. They can uncover in- 
ternal characteristics unaccessible by other means - here 
e. g. g s /g m which is not directly measurable in the most 
natural SEP layout (sample 1). At low temperature de- 
viations from the sequential behavior due to quantum 
fluctuations become clearly visible. They are described 
in close analogy to the SET by a logarithmic correction 
term of the form AG q f = a \og(ksT / (2E W )) , in qualita- 
tive agreement with the findings of Re f. I la We propose 
a reevaluation of the formulae of Ref. |l2| with parame- 
ters in accordance with our experiment so as to check the 
applicability of perturbation expansion for devices more 
complex than the SET. 
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